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Based on the nuclear attenuation data obtained by the HERMES experiment on nitrogen and kryp- 
ton nuclei, it is shown that the nuclear attenuation R\ t can be parametrised in a form of a linear 
polynomial Pi = an + ran, where r is the formation time, which depends on the energy of the 
virtual photon v and fraction of that energy z carried by the final hadron. Three widely known pa- 
rameterizations for r were used for the performed fit. The fit parameters an and a\2 do not depend 
on v and z. 
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I. INTRODUCTION 



Semi-inclusive Deep Inelastic Scattering (DIS) of lep- 
tons on nuclear targets is widely used for studies of the 
hadronization process JlJ, |2j, |3j, \4j\. It is most effective 
to observe at moderate energies of the virtual photon, 
when the formation time of the hadron is comparable 
with the nuclear radius. 

Hadronization is the process through which partons, 
created in an elementary interaction, turn into the had- 
rons that are observed experimentally. According to 
theoretical estimates the hadronization process occurs 
over length scales that vary from less than a femtome- 
ter to several tens of femtometers. Hadronization in 
a nuclear environment is interesting due to two main 
reasons. First of all, it allows to study the parameters 
governing this process on the early stage; on the other 
hand, it can provide initial conditions for the investiga- 
tion of hadronization in hot nuclear matter, which arises 
in high energy ion-ion collisions. 

The most convenient observable measured experimen- 
tally for this process is the Nuclear Attenuation (NA) 
ratio, which is a ratio of multiplicities on nucleus and 
deuterium (per nucleon). Next step is to find a variable 
which allows to present the NA ratio in the most simple 
form. 

We propose the formation time r as the best variable for 
the NA and show that the data can be parametrised in 
the form of a linear polynomial P\ = an + ra 12. Forma- 
tion time t depends on the energy of photon v and frac- 
tion of this energy z = Eh/v carried by the final hadron 
with energy Eh. Three widely known parameterizations 
for r were used for the fit procedure. The parameters 
an and a±2, obtained from the fit, do not depend on 
v and z. They are functions of the prehadron-nucleon 
and hadron-nucleon cross sections and the atomic mass 
number. 

The NA data for pions on nitrogen and for identified 
hadrons on krypton nuclei obtained by the HERMES ex- 
periment [3, 4] were used to perform the fit. 
This paper is organized as follows. Nuclear attenuation 
in absorption model is presented in the next section. In 
Section 3 it is discussed a choice of the appropriate form 



for the variable r. Section 4 presents results of the fit. 
Conclusions are given in Section 5. 



II. NUCLEAR ATTENUATION IN ABSORPTION 
MODEL 

The semi-inclusive DIS of lepton on nucleus of atomic 
mass number A is: 



li + A —> If + h + X, 



(1) 



where h(lf) are the initial (final) leptons, and h is the 
hadron observed in the final state. This process is usu- 
ally investigated in terms of the NA ratio, which is fre- 
quently defined as a function of two variables, v and z 1 



R\j{v, z) w 2(1(Ta{v, z) / Ad(j£,{v, z). 



(2) 



Usually it is investigated at precise values of one vari- 
able and average values of another. 
In case where the v - dependence is studied, 



R h M (v, (z)) = 2da A {v, {z))/Ada D {v, («)), 



(3) 



where (z) are the average values of z for each v bin. And 
for z - dependence 



R h M {(v),z) = 2da A {{v),z)/Ada D {{v),z), 



(4) 



where (v) are the average values of v for each z bin. 

In this work we adopt a model, according to which the 
origin of NA is the absorption of the prehadron (string, 
dipole) and final hadron in the nuclear medium. In that 
case NA ratio has the following form: 



R h 



M 



d Z b 



p{b,x)[W{b,x)f A - x Ux, 



(5) 



1 In fact, R/i, also depends on the photon virtuality Q 2 and on the 
square of the hadron transverse momentum pi, in respect to the vir- 
tual photon direction. However, from the experimental data, it is 



h 

M 

son with Q 2 and p| 



known that is a much sensitive function of v and z in compari- 
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where W(b, x) is the probability that neither the pre- 
hadron produced at the DIS point (b, x) (where b is im- 
pact parameter and x the longitudinal coordinate), nor 
the hadron produced at point (b, x') is absorbed by one 
nucleon in the nucleus; p is the nuclear density function 
with a normalization condition: 

J p(r)d 3 r = 1. 

For W(b, x) we use the one scale model proposed in 
Ref. H: 

W{b,x) = l-a g P q (x' - x)p{b,x')dx'- (6) 

J X 

r-OQ 

Oh I Ph{x - x)p{b,x')dx , 

J x 

where a q and ah are the inelastic cross sections for 
prehadron-nucleon and hadron-nucleon interactions, 
respectively. P q {x' — x) is the probability that on dis- 
tance x 1 — x from DIS point particle is a prehadron and 
Ph {x' — x) is the probability that particle is a hadron. 

P h (x'-x) = l-P q (x'-x). (7) 

In analogy with survival probability for a particle hav- 
ing lifetime r in a system where it travels a distance 
x' — x before decaying, P q (x' — x) can be expressed in 
form: 

P q (x'-x) =exp[-(x' -x)/t], (8) 

where r is the formation time. Substituting expressions 
for P q (x' — x) and Ph{x' — x) in eq.(6) we obtain: 



W(b, x) « 1 - a h / p(b, x')dx' + r{a h - a q )p(b, x) (9) 

J x 

w wi(b, x) + t(v, z)w 2 (b, x). 

W depends on v and z only by means of t{v, z). 
In more detail, the formation time in string models 
can be divided in two parts (see, for instance, two 
scale model presented in Ref [2, 6]). First past is the 
constituent formation time r c , which defines the time 
elapsed from the moment of the DIS untill the produc- 
tion of the first constituent of the final hadron. Second 
time interval begins with the production of the first con- 
stituent until the second one, which coincides with the 
yo-yo 2 or final hadron production. Comparison with the 
experimental data shows that in the second interval, the 
prehadron-nucleon cross section has values close to the 



hadron-nucleon cross section ah- If the difference be- 
tween these cross-sections is neglected, the model is re- 
duced to one scale model with r — r c . In case of the im- 
proved two scale model Ref [6], the prehadron-nucleon 
cross section reaches hadron-nucleon cross section value 
during a time interval r = t c + cAt, where At = zv/k, 
k is the string tension, c is the free parameter which de- 
fines from fit. In Ref [6] it is shown that c <C 1. Transi- 
tion to the one scale model takes place at c=0 and corre- 
sponds to r = t c . One should note that any complicated 
absorption string model, in some approximation, can be 
reduced to the one scale model presented in eq.(6) and 
(9). Substituting W(b, x) in R\ [ we obtain: 

R h M~fd 2 b( p(b,x)(w 1 +Tw 2 ) {A - 1) dx 

J J — oo 

w a a + ra i2 + r 2 a l3 + ■ ■ ■, (10) 

where i is the maximal power of r with which we are 
limited. The coefficients a.^ depend on A, a q , ah and nu- 
clear density. This means, that vary for different nu- 
clei. For each nucleus ay are the same for hadrons with 
equal values of a q and a h (for instance for pions 3 and 
negatively charged kaons). Although FCj^ is a polyno- 
mial of r with maximal power A — 1, it is expected that 
an > a,i2 > ai3 > ••: For fitting of the NA data we use 
three expressions for as first, second and third order 
polynomials of r : 

Pi = an + Ta 12 , (11) 
Pi = a 2 i + ra 22 + T 2 a 23 , (12) 

Ps = a 3 i + ra 32 + r 2 a 33 + T 3 a 34 . (13) 

In order to get the information on the influence of 
highest order of polynomials, also P4 polynomial was 
checked (see section 4). 

III. FORMATION TIME 

Equation (10) shows, that within our approximation, 
R'lj is a function of r only. In this section we will dis- 
cuss the physical meaning and possible expression of 
the formation time r. There are different definitions for 
the formation time. We define formation time as a time 
scale within which the prehadron-nucleon cross section 
reaches the value of the hadron-nucleon one. In the lit- 
erature there are three qualitatively different definitions 
for t. In the first extreme case it is assumed that t = 



2 The yo-yo formation means that a colorless system with valence 
contents and quantum numbers of the final hadron is formed, but 
without its "sea" partons. 



3 we do not mention the electric charge of pions, because cross sec- 
tions of differently charged pions with nucleons, which are of inter- 
est to us, are equal. 
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(Glauber approach). In the second extreme case r 3> ta, 
where ta is the nuclear radius (energy loss model 01). 
And at last, in our opinion more realistic definition of 
the formation time, as a function of v and z which can 
change from zero up to values larger than ta- Experi- 
mental data seem to confirm that for moderate values of 
v on the order of lOGeV the formation time is smaller 
than the nuclear size, i.e. the hadronization takes place 
within the nucleus. We will use the formation times ac- 
cording to the third definition mentioned above. Fol- 
lowing expressions are used: 

1. Formation time for the leading hadron [8], which fol- 
lows from the energy-momentum conservation law 



Tlead = (1 - Z)v/K, 



(14) 



where k is the string tension (string constant) with nu- 
merical value k — lGeV/fm. Indeed, eq. (14) presents 
formation time for the hadron produced on the fast end 
of the string or, which is the same, for the last hadron 
produced from string. Hadrons can be produced along 
whole length of the string. Among them, the hadrons 
produced on the fast end have better chance to avoid in- 
teractions in the nucleus. 

2. Formation time for the fast hadron, which is com- 
posed of characteristic formation time of the hadron h 
in its rest frame tq and Lorentz factor (see Ref. |5j]) 



TLo 



TO- 



Eh 



m h 



(15) 



where Eh and nih are the energy and mass of the hadron 
h, respectively. Let us briefly discuss the factor to . Some 
authors assume (see Ref. [9]) that To is a universal quan- 
tity which does not depends on the hadron type. If 
this assumtion is correct, then T Lor for the kaons is ap- 
proximately 3.5 times shorter than for the pions at the 
same values of v and z, although the experiment gives 
Km ~ R M , and it is known that ov w a K - . Such defini- 
tion seems merely strange in framework of the absorp- 
tion model. It appears more realistic that tq is propor- 
tional to nih- The reason for this is that in the string 
model the meson (baryon) is represented as a system 
consisting of a quark-antiquark (diquark) and a gluonic 
string between them. The energy of the system is trans- 
fered by the gluonic string from one parton to another 
and back. One full cycle lasts a period of nih/ k which 
we adopt as tq. Then T£ or is a universal quantity which 
does not depend on the hadron type. 
3. Formation time following from the Lund string model 
Ref. d 4 



7~Lund — 



ln(l/z 2 ) - 1 + z 2 



Z1S 
K 



(16) 



4 Note that this approximation is used only for the sake of easy read- 
ing. For numerical calculations we use the precise expression for 



TLund following from equation r Lund 
from eq.(4.21) of Ref. d. 



■ zv/k with T y taken 



One should note that all three types of formation time 
have similar behaviour with v, but different behaviour 
with z. At the values of z typical for the HERMES kine- 
matics the behavior of t defined as (14) and (16) is simi- 
lar, i.e. they are decreasing with the increase of z, while 
t defined as (15) is increasing with the increase of z. 



IV. RESULTS 

A combined fit was performed for each nucleus and 
all hadrons having equal cross sections and, conse- 
quently, identically absorbed by nuclear matter (the 
same applies to their prehadrons). It is worth remind- 
ing, that the inelastic cross sections of corresponding 
hadrons with nucleons in the moderate energy range 
{Eh ~ lOGeV) have the following values: ov = a^- = 
20mb, <jk+ = 14m6, a p — 32mb and <jp — 42mb. Com- 
bined fit was performed only for the hadrons which 
have equal cross-sections. As a result, two combined 

R h 
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FIG. 1: The values R M on nitrogen as a function of ri ea d (left two 
panels), TLor. (central two panels), TLund (right two panels). Nor- 
malized values t jr(max) for all t are used. On upper panels n + , 
on lower n~ mesons are presented respectively. Solid, dashed and 
dotted curves are results of linear, quadratic and cubic polynomial 
fits. The numerical results for linear fit are noted on upper panels. 

fits were performed for: 

1. positive and negative pions on nitrogen (26 expe- 
rimental points from [3J]); 

2. positive, negative, neutral pions and nega- 
tive kaons on krypton (63 experimental points 
from |4J]); 
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FIG. 2: 77ze values on krypton as a function ofri ea d (left four 
panels), tlot. (central four panels) , TLund (right four panels). On 
panels from upper to lower ty + , iy~ , tv° and K~ mesons, are pre- 
sented respectively. 
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FIG. 3: The values R%[ on krypton as a function of ri ea d (left three 
panels), tlot. (central three panels), r Lund (right three panels). The 
results for K + mesons are presented on upper panels, protons - on 
middle, and antiprotons - on lower. 



3. and separate fits for: 

• positive kaons on krypton (16 experimental 
points [4J]); 

• protons on krypton (16 experimental 
points (I) and 

• antiprotons on krypton (14 experimental 
points |4J]). 

The NA ratios were taken in polynomial form (11)-(13), 
and formation times (lengths) as in (14)-(16). The results 
for the reduced x 2 denoted as x 2 /d.o.f. are presented in 
Table |U One can see that for each choice of formation 
time and for each nucleus, the values of x 2 /d-o.f. are 
close for the polynomials Pi, Pi and P3, which means 
that including higher power polynomials of r will not 
essentially improve the description of data. In order 
to test this, we calculated also P4 polynomials and ob- 
tained the values of x 2 /d.o.f. close to the ones in case of 
P3. From Table U one can see that the fit gives unexpect- 
edly good values for x 2 1 d.o.f . close to unity for Ti ea( i and 
TLund, and for tl t the agreement is much worse. Expe- 
rimental points as a function of r and results of the fit 
are presented in Fig. Q] for nitrogen and in Figures |2 [3] 
for krypton. Solid points correspond to the NA ratio 
obtained from the experimental data for ^-dependence, 
open points for z-dependence. For convenience we have 
renormalized r to t / T(max), where r(max) are the max- 
imum values of r for each set of data and each choice of 
t expression. On each of the figures the linear polyno- 
mial is presented an + ra 12 with values an and a 12 = 
a\2T(max) corresponding to the best fit. Solid, dashed 
and dotted curves correspond to Pi, P2 and P3 polyno- 
mial fit. 



V. CONCLUSIONS. 

Based on the performed studies one can assume 
that NA ratio is a function of formation time r 
only. 

Three expressions for NA ratio as a first, second 
and third order polynomials of t were used for 
these studies, and the results do not show remark- 
able difference in their ability to describe the HER- 
MES data. 

The performed analysis with three different ex- 
pressions for the formation time r shows that 
based on the obtained values of the reduced x 2 
(see Table 1), as well as on visual shape of the 
curves in Figures [TJ|2] and |3) one can conclude that 
the two of them - Ti ea d and TLund are quite appro- 
priate, and the expression for the formation time 
in form of tlot is ruled out. 

The performed combined fits are based on the 
HERMES NA data using: tt + , it~ for nitrogen; n + , 
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7r , 7T , K for krypton; separate fits for K + , pro- 
tons and antiprotons for krypton. 

• Because of the fact that the linear function gives 
about the same values x 2 /d.o.f. values as second 



and third order one, we assume that in our ap- 
proach the parameterization by linear function is 
well enough. 
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TABLE I: x 2 /d.o.f. from polynomial fit. Details see in text 
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